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Abstract 

The Boltzmann equation for inelastic Maxwell models is used to analyze non- 
linear transport in a granular binary mixture in the steady simple shear flow. 
Two different transport processes are studied. First, the rheological proper- 
ties (shear and normal stresses) are obtained by solving exactly the velocity 
moment equations. Second, the diffusion tensor of impurities immersed in 
a sheared inelastic Maxwell gas is explicitly determined from a perturbation 
solution through first order in the concentration gradient. The correspond- 
ing reference state of this expansion corresponds to the solution derived in 
the (pure) shear flow problem. All these transport coefficients are given in 
terms of the restitution coefficients and the parameters of the mixture (ratios 
of masses, concentration, and sizes). The results are compared with those 
obtained analytically for inelastic hard spheres in the first Sonine approxima- 
tion and by means of Monte Carlo simulations. The comparison between the 
results obtained for both interaction models shows a good agreement over a 
wide range values of the parameter space. 
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I. INTRODUCTION 



Granular media under rapid flow conditions are usually modeled by a fluid of hard spheres 
with inelastic collisions. In the simplest version, the grains are taken to be smooth so that 
the inelasticity is only characterized through a constant coefficient of normal restitution. 
For a low-density gas, the Boltzmann equation has been conveniently generalized 1 to ac- 
count for the inelasticity of binary collisions and the Navier-Stokes transport coefficients 
have been obtained in terms of the restitution coefficient. 2-5 As in the case of elastic col- 
lisions, these transport coefficients verify a set of coupled linear integral equations which 
are solved approximately by using the leading terms in a Sonine polynomial expansion. 
The corresponding Sonine results are in general in good agreement with selected tests using 
molecular dynamics 6 and Monte Carlo simulations. 4 ' 7,8 However, beyond the Navier-Stokes 
regime (small spatial hydrodynamic gradients), it is quite intricate to get explicit results for 
inelastic hard spheres (IHS) and one has to resort to alternative approaches. 

One of the main mathematical difficulties in solving the Boltzmann equation for hard 
spheres (even in the elastic case) comes from the form of the collision rate, which is pro- 
portional to the magnitude of the relative velocity of the two colliding spheres. In the case 
of elastic fluids, a possible way to overcome this problem is to assume that the particles 
interact via the repulsive Maxwell potential (inversely proportional to the fourth power of 
the distance). For this interaction model, the collision rate is independent of the relative 
velocity and this allows for a number of nice mathematical properties of the Boltzmann 
collision operator. 9 Thanks to this simplification, nonlinear transport properties can be ex- 
actly obtained 10 from the Boltzmann equation for Maxwell elastic molecules and, when 
properly reduced, they exhibit a good agreement with results obtained for other interaction 
models. 11 In the context of inelastic gases, the Boltzmann equation for inelastic Maxwell 
models (IMM) was also introduced about two years ago. 12-14 The IMM share with elastic 
Maxwell molecules the property that the collision rate is velocity independent but their 
collision rules are the same as for IHS. As a consequence, these IMM's do not describe real 
particles since they do not interact according to a given potential law. For this reason the 
model is usually referred to as pseudo-Maxwellian model. 13 Nevertheless, the IMM keeps the 
qualitatively correct structure and properties of the nonlinear macroscopic equations and 
obey Haff's law. 15 

The Boltzmann equation for IMM has received a great attention in the last few years, 
especially in the study of overpopulated high energy tails in homogeneous states. 15-18 The 
existence of high energy tails of the Boltzmann equation is common for IHS and IMM, 
although this general qualitative agreement fails at a quantitative level. For inhomogeneous 
states, the Chapman-Enskog method has been recently applied to the Boltzmann equation to 
get the Navier-Stokes transport coefficients. 19,20 The comparison with the results derived for 
the IHS model 2 shows in general significant discrepancies in the dependence of the transport 
coefficients on the restitution coefficient. The IMM has been also used to get the shearing 
stress of a granular material under shear flow. 21 All the above results refer to monocomponent 
gases. Much less is known for inelastic Maxwell mixtures. In the case of multicomponent 
gases, Marconi and Puglisi 22 have analyzed energy nonequipartition in the free cooling and 
driven states for a one-dimensional system, while Ben-Naim and Krapivsky 23 have studied 
velocity statistics of an impurity inmersed in a uniform granular fluid. However, to the best 
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of my knowledge, no previous study on transport properties in inelastic Maxwell mixtures 
has been carried out. 

The aim of this paper is to analyze nonlinear transport in a binary mixture in the frame- 
work of the Boltzmann equation for IMM. Two different transport processes will be studied. 
First, the rheological properties (shear stress and normal stress differences) of a granular 
binary mixture subjected to the simple shear flow will be explicitly determined. Second, the 
elements of the tracer diffusion tensor characterizing the diffusion of impurities in a granular 
gas under shear flow will be obtained in the first order of the concentration gradient. In 
both situations, given the coupling between dissipation and the shear rate, the mixture is 
far away from equilibrium and the Navier-Stokes constitutive equations 2 ' 19 ' 20 do not apply. 
This makes the analytical treatment of both nonequilibrium situations difficult. An addi- 
tional complication arises from the large number of parameters governing the dynamics of 
the system, including three independent restitution coefficients. However, the tractability of 
the Boltzmann collision integrals for IMM allows one to determine all the above nonlinear 
transport properties by solving exactly the hierarchy of moment equations. The results de- 
rived for IMM are then compared with those recently obtained for IHS in the leading Sonine 
approximation 24 ' 25 and by means of Monte Carlo simulations. 24,26 Since the strength of the 
shear rate is arbitrary the test of the capability of IMM to capture the relevant behavior of 
IHS is quite stringent. The comparison performed here shows that in general the agreement 
between both interaction models is quite good, even for strong dissipation. 

The plan of the paper is as follows. In Sec. II the Boltzmann equation for IMM and the 
macroscopic conservation laws are introduced. Section III deals with the simple shear flow 
problem for inelastic Maxwell mixtures. The elements of the pressure tensor are exactly 
obtained in terms of the restitution coefficients and the parameters of the mixture (masses, 
sizes, and concentrations). In particular, the results show that both species have different 
temperatures so that the energy is not equally distributed between both species (breakdown 
of energy equipartition). The diffusion tensor of impurities in a Maxwellian sheared fluid is 
obtained in Sec. IV. The paper is closed in Sec. V with some concluding remarks. 



II. THE BOLTZMANN EQUATION AND THE INELASTIC MAXWELL MODEL 

Let us consider a binary mixture of inelastic Maxwell gases at low density. In the simplest 
version, the Boltzmann equation for IMM 12 ' 15 ' 18,19,23 can be obtained from the Boltzmann 
equation for IHS by replacing the rate for collisions between particles of species r and s 
by an average velocity-independent collision rate, which is proportional to the square root 
of the temperature. This means that a random pair of colliding particles undergo inelastic 
collisions with a random impact direction. With this simplification, the set of nonlinear 
Boltzmann kinetic equations become 



{Ft + Vl ' V ) /r(r ' Vi; t] = ^ Jrs [vi|/ ^' fs ^ ] ' 



(1) 



where / r (r, Viji) is the one-particle distribution function of species r (r = 1,2) and the 
Boltzmann collision operator J rs [vi|/ r , f s ] describing the scattering of pairs of particles is 

Jrs[vi\f r ,fs\ = -^77- / dv 2 / da \a~sfr(r,v[,t)f s (r,v' 2 ,t) - f r (r, Vi, t)f s (r, v 2 , *)1 • (2) 

U s l L f : J J L J 
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Here, n s is the number density of species s, w rs is an effective collision frequency (to be 
chosen later) for collisions of type r — s, Qd — 27r d / 2 /Y{d/2) is the total solid angle in d 
dimensions, and a rs < 1 refers to the constant restitution coefficient for collisions between 
particles of species r with s. In addition, the primes on the velocities denote the initial 
values {v' 1? v 2 } that lead to {vi, v 2 } following a binary collision: 

v i = Vi - n„ (l + a^ 1 ) (ct • gi 2 )3", v 2 = v 2 + p rs (l + a" 1 ) (5- • g 12 )5- , (3) 

where g 12 = Vi — v 2 is the relative velocity of the colliding pair, a is a unit vector directed 
along the centers of the two colliding spheres, and 

fl rs = ■ . 4 

m r + m s 

There is another more refined version of the inelastic Maxwell model 13 ' 14 ' 27 where the collision 
rate has the same dependence on the scalar product (<r • gi 2 ) as in the case of hard spheres. 
The corresponding Boltzmann equation can be proved to be equivalent to Eq. (1) except 
that the collision laws (3) must be changed. However, both versions of IMM lead to similar 
results in problems as delicate as the high energy tails. 15 Therefore, for the sake of simplicity, 
here I will consider the version given by the collision rules (3). As will be shown later, this 
model compares quite well with the results obtained from the IHS model. 

At a hydrodynamic level, the relevant quantities in a binary mixture are the number 
densities n r , the flow velocity u, and the "granular" temperature T. They are defined in 
terms of moments of the distribution f r as 



71, 



= J dv/ r (v), pu = J2pr- U r = J2j dvm r vf r (v), (5) 



nT = 5>T r = £/ dv^V 2 f r (v), (6) 

where p r = m r n r , n = n x + n 2 is the total number density, p = p\ + p 2 is the total mass 
density, and V = v — u is the peculiar velocity. Equations (5) and (6) also define the flow 
velocity u r and the partial temperature T r of species r, the latter measuring the mean kinetic 
energy of species r. 

The collision frequencies w rs can be seen as free parameters in the model. Its dependence 
on the restitution coefficients a rs can be chosen to optimize the agreement with the results 
obtained from the Boltzmann equation for IHS. Of course, the choice is not unique and may 
depend on the property of interest. 

The collision operators conserve the particle number of each species and the total mo- 
mentum, but the total energy is not conserved. This implies that 

£ / dv\m r V 2 J rs [v\fr, fs] = -fnTC , (7) 

r,s J 

where ( is identified as the "cooling rate" due to inelastic collisions among all species. At a 
kinetic level, it is also convenient to discuss energy transfer in terms of the "cooling rates" 
£ r for the partial temperatures T r . They are defined as 
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Cr — ^2(rs, (8) 

s 

where 

(rs = --r-TF I dvm r V 2 J rs [v\f r ,f s ] . (9) 



dn r T 7 

The total cooling rate ( can be expressed in terms of the partial cooling rates ( r as 

C = T- 1 Y J XrT r C, r , (10) 

r 

where x r = n r /n is the mole fraction of species r. 

The macroscopic balance equations follow from the Boltzmann equations (1) by taking 
velocity moments. They are given by 

^7 * 

D t n r + n r W -u + — — = 0, (11) 

m r 

Au + p- 1 VP = 0, (12) 

£> t T - - £ — + ^- (V • q + P : Vu) = -CT . (13) 
n r m r an 

In the above equations, D t = d t + u ■ V is the material derivative, 

j r = m r fdvVf r (v) (14) 
is the mass flux for species r relative to the local flow, 

P = E / dvm r VV/ r (v) (15) 
is the total pressure tensor, and 

q = E / ^v|m r \/ 2 V/ r (v) (16) 

is the total heat flux. 

The main advantage of using IMM is that a velocity moment of order k of the Boltzmann 
collision operator only involves moments of order less than or equal to k. 10 ' 27 This allows 
one to determine the Boltzmann collision moments without the explicit knowledge of the 
velocity distribution function. In general, beyond the linear hydrodynamic regime (Navier- 
Stokes order), the above property is not sufficient to exactly solve the hierarchy of moment 
equations due to the free-streaming term of the Boltzmann equation. Nevertheless, there 
exist some particular situations (such as the simple shear flow problem) for which the above 
hierarchy can be recursively solved. The first few moments of the Boltzmann collision 
operator J rs [f r , fs] have been evaluated in Appendix A. They can be written as 
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/ dvm r VJ rs [f r , f s ] = -^/i sr (l + a rs ) (pj r - p r j s ) , 



(17) 



/w 
dvm r WJ rs [f r , f s ] = ^/i sr (l + a rs ) {2p s P r - (j r j s + j j r 



d + 2 
d 



H sr (l + a rs ) [p s P r + p r P s - (iris + j s j r ) 



+ 



(PrPs + PsPr) - ir is 
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where 



= J dvm r W/ r , 



(18) 



(19) 



and p r = n r T r = trP r /d is the partial pressure of species r. The relationship between the 
collision frequency w rs and the cooling rate ( rs can be obtained by taking the trace in Eq. 
(18): 



_ 

Srs ^ P>sr\i- ~r OVsJ 



1 - -7r(l + «r S ) T 1 "J Jr • Js 

2 S dpsPr 



where 



^217 

and 7 = Ti/T 2 is the temperature ratio. 



e ^ = 1 + ^(7-1) ^ = i + a;i ( 7 -i) 



(20) 



(21) 



III. RHEOLOGICAL PROPERTIES UNDER SIMPLE SHEAR FLOW 



As said in the Introduction, one of the main objectives of the present paper is to evaluate 
the rheological properties of an inelastic Maxwell binary mixture subjected to the simple 
shear flow. At a macroscopic level, this state is characterized by a constant linear velocity 
profile u = ui = u 2 = a • r, where the elements of the tensor a are a^- = aSi X 5j y , a being 
the constant shear rate. In addition, the partial densities rii and the granular temperature 
T are uniform, while the mass and heat fluxes vanish by symmetry reasons (j r = q = 0). 
As a consequence, the (uniform) pressure tensor P is the only nonzero flux in the problem. 
On the other hand, according to the energy balance equation (13), the temporal variation 
of the granular temperature arises from the balance of two opposite effects: viscous heating 
and collisional cooling. When both mechanisms cancel each other, a steady state is achieved 
and the temperature remains constant. In that case, the shear stress P xy and the cooling 
rate ( are related by 



aR 



xy 



(22) 



where p = trP/d = nT is the pressure. Our goal here is to determine the pressure tensor P 
in this steady state. 
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The simple shear flow becomes spatially uniform when one refers the velocities of the 
particles to a frame moving with the flow velocity u: f s (r, v) — > / S (V). Consequently, the 
corresponding steady Boltzmann equations (1) read 



d 

- a ij V 3gy,fl = Ju[fl, h\ + ^12 [/l, H 



(23) 



d 

-aijVj—h = J 22 [f 2 , f 2 ] + J21U2, fi] ■ 



(24) 



The total pressure tensor is P = Pi + P2, where the partial pressure tensors P r (r = 1, 2) 
are defined by Eq. (19). The elements of these tensors can be obtained by multiplying the 
Boltzmann equations (23) and (24) by m r VV and integrating over V. The result is 

a ik Pi,kj + a jk Pi,ki + B u P 1Aj + B 12 P 2tij = (Aupx + A u p 2 ) (25) 

where use has been made of Eq. (18) (with j r = 0) and I have introduced the coefficients 



A 



M = 2(^) (1 + ail) + ^ 2l(1 + ai2) ' 



(26) 



4 W 12 Pi 2 /-, , x2 

Al2 = ZT^^ 2l(1 + ai2) ' 



(27) 



= „T n ^ (l + «n)( rf + 1 " «n) + J? n «s toiO- + "12) [d + 2- /i 2 i(l + «i 2 )] , (28) 



d(d + 2) 



d(d+2)' 



B 



12 



2 . 

-1^12. 
a 



(29) 



A similar equation can be obtained for P 2 , by just making the changes 1 2. From Eq. 
(25), in particular, one gets 



aP 1 



xy 



d 



(30) 



In order to have explicit expressions one still needs to fix the parameters w rs . The most 
natural choice to find good agreement with the IHS results is to adjust the cooling rates 
( rs for IMM, Eq. (20), to be the same as those obtained for IHS. Given that the cooling 
rates for IHS are not exaclty known in the simple shear flow problem, I will take for ( rs the 
cooling rate of IHS evaluated at the first Sonine approximation: 24 



IHS 

rs 



2tt 



y/7ld 



n s /i sr o- sr v 



6 r + 6 S 
9 r 9. q 



1/2 



(1 + a r 



1 Psr /-. . \9 r ^ S 

1 - — (1 + a rs )— - — 



(31) 



where a rs = (a r + <r s )/2 and v = y / 2T(m i + m 2 )/mim 2 is a thermal velocity defined in 
terns of the granular temperature of the mixture T. Therefore, according to Eq. (20), the 
collision frequencies w rs of IMM are given by 
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w r . = ^n.o£\\^£.) , (32) 

where use has been made of the fact that j r = in the simple shear flow state. Henceforth, 
I will present the results of IMM with this choice. 

Equation (25) along with its counterpart for P 2 constitute a linear system of equations for 
the total pressure tensor P. In order to express the solution of this system, it is convenient 
to introduce dimensionless quantities. As done in a previous study for IHS, 24 we reduce the 
shear rate a with respect to an effective collision frequency 

ild na d l2 l v . (33) 



2v^F 

Thus, we introduce the reduced cooling rates (* = (, r jv, the reduced shear rate a* = a/u, 
and the reduced partial pressure tensors P* = P r /x r p. The reduced total pressure tensor is 
P* = P/p = x x P\ + x 2 P* 2 . It is worthwhile remarking that, for given values of the parameters 
of the mixture, a* and P* are only functions of the restitution coefficients an, a 22 , and a 12 . 

In reduced units, the solution to Eq. (25) gives the pressure tensor P* in terms of a*, 
the temperature ratio 7, the restitution coefficients and the parameters of the mixture. The 
dependence of a* on a rs can be obtained from the energy equation (22) as 

a * = = d gffiCj + ^272(2* ^ 

2 P* y 2 XiP* xy + X 2 P 2xy 

where I have introduced the temperature ratios 71 = T\jT = (/i2i#i) 1 and 72 = T 2 /T = 
(A*i2^2) _1 - When the dependence of P* and a* on the temperature ratio 7 is known, the 
latter can be obtained by solving Eq. (30) (or its counterpart for the species 2) 

a* p* 

7 = J^F- (35) 

SI *2,xy 

Once Eq. (35) is solved, one gets the rheological properties (as measured by the elements 
of the pressure tensor) of the mixture in terms of the parameters of the problem. In the 
elastic limit (an = a 22 = a 22 = 1), the cooling rate ( = and the temperature increases 
in time due to viscous heating (aP xy 7^ 0). Consequently, for long times, the reduced shear 
rate a*(t) = a/u(t) tends to zero and one recovers the equilibrium results, P*j = 5ij. Beyond 
the elastic limit, the reduced pressure tensor has a complex dependence on the restitution 
coefficients and the parameters of the mixture. Before considering the general case, let us 
study separately some interesting limit cases. 



A. Mechanically equivalent particles 

This simple case corresponds to m x = m 2 , <j\ = <j 2 , and an = a 2 2 = a i2 = a. In this 
limit, Eqs. (25), (34), and (35) yield 7 = 1, P* = P* = P*, and the nonzero elements of the 
(reduced) pressure tensor can be written as 



P*yy -PI-'-'- P* dd ~ 2 d+ + x _ a i ( 36 ) 
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xy 



d 2 (d + 2)(l-a 2 ) 
~\ 8(d+l-a) 2 ' 



(37) 



PL = d - {d - i) P * yy , (38) 

Expressions (36)-(38) are independent of the choice of the collision frequency w rs = w. 

Recently, Cercignani 21 has analyzed the simple shear flow problem of a monocomponent 
granular gas considering the more refined version of the Maxwell model indicated in Sec. 
II. The results obtained by Cercignani 21 in the (steady) shear flow problem for a three 
dimensional system leads to the absence of the normal stress differences (P* x = P* y = P* z = 
1) and the shear stress is given by 

p * = -f^- < 39 > 

The elements of the pressure tensor of IHS have been explicitly obtained from the Boltz- 
mann equation in the first Sonine approximation. 24 ' 25 For the sake of completeness, their 
expressions are displayed in Appendix B [cf. Eqs. (B5-(B7)]. It is apparent that these ex- 
pressions differ from those found here, Eqs. (36)-(38). However, for practical purposes, the 
discrepancies between both interaction models are quite small, even for moderate values of 
a. As an ilustration, Fig. 1 shows a comparison of the results obtained for IMM and IHS 
for the (reduced) elements of the pressure tensor for d = 3 . The results (symbols) obtained 
from Monte Carlo simulations for IHS 24 as well as the predictions for the shear stress given 
by Eq. (39) are also included. We observe that the discrepancies between IHS and IMM 
are very small, although they slightly increase as the restitution coefficient decreases. In 
addition, the agreement of both interaction models with simulation data is quite good. On 
the other hand, the results derived in ref. 21 exhibit significant discrepancies with Monte 
Carlo results, except in the low-dissipation region. 



B. Tracer limit 

When one of the species, say species 1, has a vanishing mole fraction, it acts as a tracer 
species. In this case, the state of the excess component 2 is not affected by the presence of 
the tracer particles and so the expression of is the same as the one obtained in the single 
Eqs. (36)-(38), with a = 0:22- The corresponding expressions for the elements of 
the tracer pressure tensor can be obtained from Eq. (25) by taking carefully the limit 
xi — > 0. After some algebra, one gets 

TT _|_ ff P* 

p* _ p * _ _ p* _ _ z,yy (An\ 
r i,yy — r i,zz — ■ ■ ■ — r ldd — — , {^U) 

a *p* _ up* 
p* = » l,yy nr %. 
l,xy q 

Pi,xx = *y-{d- 1)PT >W , (42) 
9 



(41) 



where 



(43) 



1/2 

./7777m... f i^il 

d{d + 2) 




(l + M ^2i(l + ai 2 ) -d-2], (44) 




H= W^2)^{—) + (45) 
and the reduced shear rate is 

.2 _ (°2\ 2(d - 1} 2(a'+l-a 22 ) 2 (l-o| 2 ) 

~w d 2 (d + 2) • l4bj 

Here, = miT 2 /m 2 Ti is the mean square velocity of the gas particles relative to that of 
the tracer particles. The temperature ratio 7 is obtained from Eq. (35) where the (reduced) 
cooling rates Q an d Ci i n the tracer limit are given by 

^ = — {—) (l-« 22 ), (47) 



Cl = -1^21 —5— (1 + «12) 



d' 



(48) 



C. General case: comparison with IHS results 

Here, I compare the results obtained for IMM with those found for IHS by using the lead- 
ing Sonine approximation 24 as well as by performing Monte Carlo simulations. 24 ' 26 Specifi- 
cally, I want to explore the dependence of the temperature ratio 7 and the nonzero elements 
of P* on the restitution coefficients a rs , the mass ratio mi/m 2 , the concentration ratio 7ii/n 2 , 
and the ratio of sizes 0"i/er 2 . Given the high dimensionality of the parameter space, for the 
sake of concreteness, henceforth I will assume that the disks or spheres are made of the same 
material, i.e., an = a 22 = a i2 = a. This reduces the parameter space to four quantities. 

The temperature ratio 7 = T\/T 2 measures the breakdown of energy equipartition. Due 
to the coupling between dissipation and the shear rate in the (steady) simple shear flow 
state [cf. Eq. (34)], the mixture departures from the Navier-Stokes regime as the restitution 
coefficient a decreases. Thus, energy nonequipartition is expected except when a — > 1. 
The coexistence of two granular temperatures has been recently confirmed by molecular 
dynamics simulations carried out in sheared granular mixtures of hard disks 28 and by some 
experiments 29 of vibrated granular mixtures. Related findings have been reported by some 
authors by using kinetic theory tools in the freely cooling state for IHS 30 ' 31 and IMM. 22,23 
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As expected, the results obtained here show that in general the kinetic temperatures of 
the mixture are different (7 ^ 1). The dependence of 7 on the parameters of the problem 
is very complex. In Fig. 2 the temperature ratio T\jT 2 is plotted versus the restitution 
coefficient a for an equimolar mixture (ni = n 2 ) of hard spheres (d = 3) with a±/a 2 = 1 
and two different values of the mass ratio: mi/m 2 = 2 and 10. We clearly see that the 
predictions for IMM agree quite well with the ones obtained for IHS in the first Sonine 
approximation and that the agreement of both approaches with simulation is excellent, even 
for moderate dissipation. It is also apparent that the extent of the equipartition violation 
is greater when the mass disparity is large. To illustrate this point, in Fig. 3 we consider 
an equimolar mixture of hard disks (d = 2) of constant density, i.e., mi/m 2 = {p\j(y 2 ) 1 . 
Three different values of the restitution coefficient have been studied: a = 0.95, 0.9, and 0.8. 
Again, the IMM results reproduce quite well the IHS predictions (theory and simulation) 
over the entire range of values of size and mass ratios considered. Further, the behavior of 
Ti/T 2 obtained here for dilute systems is qualitatively similar to that found in molecular 
dynamics simulations 28 for finite-density systems. Thus, for instance, at a given value of a, 
the granular energy of the larger particle (say for instance, species 1) increases relative to 
that of the smaller particle as the ratio o\jo 2 (or equivalently, m\jm 2 ) increases. 

The dependence of the nonzero elements of the pressure tensor P* on the restitution 
coefficient a is illustrated in Figs. 4 and 5 for two different mixtures. In Fig. 4, we consider 
a mixture of hard spheres (d = 3) with n\jn 2 = oxja 2 = 1 and m 1 /m 2 = 2, while Fig. 5 
corresponds to a mixture of hard disks (d = 2) with n\jn 2 = 1, m\jm 2 = ((X1/02) 2 = 10. We 
see that the agreement between the results derived for IMM and IHS is again in general quite 
good, especially in the case of the shear stress P* y which is the most relevant rheological 
property in a shear flow problem. In fact, in the two-dimensional case, the exact results 
obtained here for IMM are closer to the simulation results than those obtained by using the 
first Sonine approximation for IHS. 24 It must be remarked that both theories only predict 
normal stress differences in the plane of shear flow (P* x > P* y = P* z = ■ ■ ■ = P^ d ), while 
the simulation results also shows anisotropy in the plane orthogonal to the flow velocity 
(P* z > Pyy)- However, these relative normal stress differences are very small and tend to 
zero as the restitution coefficient decreases. 

IV. DIFFUSION OF IMPURITIES UNDER GRANULAR SHEAR FLOW 

The second problem analyzed in this paper refers to the diffusion of impurities in a 
granular Maxwell gas subjected to the simple shear flow. This problem has also been recently 
analyzed by the author 25 in the framework of IHS. In the tracer limit (x\ — > 0), the state 
of the excess component is not affected by the presence of the tracer particles so that its 
velocity distribution function f 2 verifies a closed Boltzmann equation. As a consequence, 
the pressure tensor P 2 of the granular gas is given by Eqs. (36)-(38). Further, given that 
the mole fraction of the tracer particles x\ is very small, one can neglect their mutual 
interactions in the kinetic equation for j\ and assume that f\ verifies a (linear) Boltzmann- 
Lorentz equation. The diffusion process is induced in the system by the presence of a weak 
concentration gradient Vx\. Given that the strength of the shear rate a is arbitrary, the 
mass flux ji (which is generated by the gradient Vxi) can be modified by the presence of the 
shear flow. Under these conditions, a diffusion tensor D is required to describe the diffusion 
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process instead of a single diffusion coefficient. The determination of the nonzero elements 
of this tensor is the goal of this Section. 

In the local Lagrangian frame, the kinetic equation governing the evolution of the velocity 
distribution function f 1 reads 

pi ~ ^y^h + (V, + ayr^/i = Jn[V\fi, / 2 ] , (49) 

where here the derivative djdri is taken at constant V. Assuming that the mole fraction 
x\ is slightly nonuniform, we solve Eq. (49) by means of a perturbation expansion around a 
nonequilibrium state with arbitrary shear rate. Thus, we write 

h = ff + A (1) + ■ ■ ■ , (50) 

where f[® is of order k in V#i but applies for arbitrary degree of dissipation since this 
distribution retains all the orders in a*. The solution (50) qualifies as a normal solution since 
all the space and time dependence of f\ occurs entirely through #i(r; t) and their gradients. 
The zeroth-order approximation f[ ^ corresponds to the simple shear flow distribution but 
taking into account now the local dependence on the mole fraction x±. Although the explicit 
form of f[ ^ is not exactly known, only the knowledge of its second-degree moments is 
necessary to get the diffusion tensor. The explicit expressions of the (reduced) pressure 
tensor P\ are given by Eqs. (40)-(45). 

The kinetic equation for can be obtained from the Boltzmann-Lorentz equation (49) 
by collecting all the terms of first order in Vxi: 

|/i (0) - a ^^A (1) + {Vi + ay^yf = JuWlfP, h] ■ (51) 

The determination of the diffusion tensor from Eq. (51) follows similar mathematical steps 
as those made in the case of IHS. I refer the reader to ref. 25 for more specific details and 
here only some intermediate results will be displayed. Using the balance equation for the 
mass density of impurities, Eq. (51) becomes 

(^4 + a)/! 1 '^(V.V, 1 ), (52) 

where A is the Boltzmann-Lorentz collision operator 

A/i (1) = Ji 2 [V|A (1) ,/ 2 ]. (53) 



The mass flux is defined as 



The components of can be obtained by multiplying Eq. (52) by m{V and integrating 
over V. The result can be written as 

(a + XH)-^ = -p 2 Pl-Vx u (55) 



jS 1 )= mi |rfvVA (1) . (54) 
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where use has been made of Eq. (17) (with j 2 = 0) and II is the d x d unit tensor. The 
coefficient A is 



A = — /i 2 i(l + ai 2 ) 

o (d-l)/2 /i i rt\ V 2 

= 5 rxw n2/i2l0r " 2 " 1(2T2/m2)1/2(1 + au) (~H ' (56) 

Note that, according to Eq. (20), wu is still given by Eq. (32) since j 2 = 0. The solution to 
Eq. (55) can be written in the form 

3? = -D-V Xl , (57) 
with the elements of the diffusion tensor D given by 

D ij = P i (8* ~ K k , (58) 



A \ lk A J W 

Equation (58) provides an explicit expression of the tracer diffusion tensor of a granular 
Maxwell mixture under shear flow. In the absence of shear rate (a = 0), as well as in the 
elastic limit ( a 22 = «i2 = 1), one has = D 5ij, where 



d T(d/2) m 1 (m 1 +m 2 ) 
D °- Vlvr^V^VV m 2 T2 (59) 

is the tracer diffusion coefficient of a gas of Maxwell molecules. 32 As the restitution coefficient 
decreases, rheological effects become important and the elements of the diffusion tensor are 
different from the one obtained in the equilibrium case. The dependence of the diffusion 
coefficients on the restitution coefficients a 22 and a 12 as well as on the mass ratio m\jm 2 and 
the size ratio 0\ja 2 is highly nonlinear. As happens for elastic fluids, 33 Eq. (58) shows that 
diffusion under simple shear flow is a very complex problem due basically to the anisotropy 
induced in the system by the shear flow. 

To illustrate the dependence of the elements of the diffusion tensor on the parameters of 
the problem, let us consider the three dimensional case id — 3). According to Eq. (58), D xz = 
D zx = D yz = D zy = 0, in agreement with the symmetry of the problem. Consequently, there 
are five relevant elements: the three diagonal and two (D xy , D yx ) off-diagonal elements. In 
general, D xx ^ D yy = D zz and D xy ^ D yx . The off-diagonal elements are negative and 
measure cross effects in the diffusion of impurities in a sheared granular gas. The results 
obtained 25 for IHS [cf. Eqs. (B8)-(B13) of Appendix B] also predict anisotropy in the plane 
perpendicular to the shear flow (D yy ^ D zz ), although this anisotropy is much smaller than 
the one observed in the plane of shear flow (as measured by the difference D xx — D yy ). In Fig. 
6, the self-diffusion problem is considered, i. e., when the tracer particles are mechanically 
equivalent to the gas particles. This figure shows D* x -D* y , (D* xx +D yy +D zz ) /3 = (l/3)D* kk , 
D* y , and D* x as functions of the restitution coefficient a = a i2 = a 22 . Here, D*j = D^/Do, 
with D given by Eq. (59). We see that the deviation from the functional form for elastic 
collisions is quite important even for moderate dissipation. The figure also shows that 
the anisotropy of the system, as measured by the difference D* xx — D yy grows with the 
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inelasticity. It is worthwhile remarking that the general qualitative behavior of the self- 
diffusion tensor on dissipation agrees well with molecular dynamics simulations carried out 
by Campbell, 34 at least for the lowest solid fraction considered in his simulations. With 
respect to the comparison with the results for IHS, 25 we see again good agreement although 
the discrepancies between both interaction models slightly increases as a decreases. In fact, 
these discrepancies are more significant than those found for the rheological properties in 
the simple shear flow problem. Finally, Fig. 7 shows the dependence of D*j on a in the 
case mi/m 2 = 2 and <7i/<7 2 = 2 for IMM and IHS. Similar conclusions to those made in the 
self-diffusion case are obtained here for mixtures of species with different (but not disparate) 
mass and size. It must be pointed out that the discrepancies for the diffusion tensor between 
both interaction models turn out to be significant as the disparity of masses or sizes increases. 

V. CONCLUDING REMARKS 

Needless to say, the analysis of nonlinear transport phenomena (i.e., when the hydrody- 
namic gradients are not necessarilly small ) from the Boltzmann equation for inelastic hard 
spheres (IHS) is a formidable task. For this reason, to get explicit results one usually consid- 
ers the leading order in a Sonine polynomial expansion of the velocity distribution function. 
All the above difficulties increase when one studies polydisperse systems (granular mixtures) 
since not only is the number of transport coefficients larger than for a monomponent gas 
but they are also functions of more parameters such as the concentrations, masses, sizes, 
and the coefficients of restitution. A possible way to partially overcome these problems is to 
consider interaction models (such as the inelastic Maxwell model) that simplify the math- 
ematical structure of the Boltzmann collision operator. For elastic fluids, this strategy has 
been shown to be very fruitful since the exact results obtained for Maxwell molecules for 
many interesting transport properties (both linear and nonlinear) agrees quite well with hard 
spheres and other interaction potentials when they are convenienlty nondimensionalized. 10 

In the Boltzmann equation for inelastic Maxwell models (IMM), the collision rate of 
IHS is replaced by an effective collision rate independent of the two colliding particles. The 
simplicity of this interaction model allows us for instance, to get exactly the moments of 
the Boltzmann collision operator without the explicit knowledge of the velocity distribution 
function. The cost of sacrificing physical realism can be compensated by the amount of exact 
analytical results obtained from IMM. However, in order to assess the degree of usefulness 
of IMM as a prototype model for the description of granular media, one needs to compare 
the predictions based on IMM with those obtained in the framework of the IHS model. 
In this paper I have carried out such a comparison at the level of nonlinear transport in 
two different but related nonequilibrium problems. First, I have explicitly determined the 
rheological properties (pressure tensor) of a binary inelastic Maxwell mixture under simple 
shear flow. Then, by taking the above solution as the reference state, I have obtained the 
diffusion tensor of impurities immersed in a sheared inelastic Maxwell gas. In the context 
of IMM, both expressions are exact in the sense that they apply for arbitrary values of the 
restitution coefficients and the parameters characterizing the mixture (masses, sizes, and 
concentrations). 

Previous studies of such nonlinear transport coefficients have been carried out 24,25 by 
the author in the case of IHS. As said above, these results 24 ' 25 are approximate since they 
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are obtained by retaining the first terms in a Sonine polynomial expansion of the velocity 
distribution functions. However, the accuracy of these results has been confirmed by Monte 
Carlo simulations 24,26 of the Boltzmann equation. To compare the results derived from IMM 
and IHS, the collision frequencies w rs appearing in the Boltzmann equation for IMM [see Eq. 
(2)] need to be fixed. Here, to make contact with IHS, I have chosen w rs to reproduce the 
cooling rates ( rs of IHS. With this choice, the results derived here for IMM for the pressure 
tensor and the diffusion tensor compare quite well with the ones previously derived for IHS 
over a wide range of values of dissipation and parameters of the mixture. Given the high 
dimensionality of the parameter space explored, the comparison carried out in this paper 
can be considered as a good test to gauge the reliability of IMM to reproduce the main 
trends observed for IHS in the context of granular mixtures. 

Finally, it must be noted that the good agreement found here contrasts with the results 
recently derived for the Navier-Stokes (NS) transport coefficients, 19 where it was shown 
that the dependence of the transport coefficients of IMM on the restitution coefficient only 
compares well at a qualitatively level with the one obtained for IHS. However, the situation 
analyzed in ref. 19 is completely different to the one studied here since the Navier-Stokes 
coefficients were obtained in the first order of the Chapman-Enskog expansion around a 
time-dependent state (homogeneous cooling state). In this sense, it seems that the results 
found in this paper for IMM may be sufficiently representative of the trends observed for IHS 
in some particular situations (such as the simple shear flow problem), and this agreement 
must be taken with caution when one explores other situations. This fact stimulates the 
search for exact solutions for IMM which can be confronted with the results obtained for 
IHS by using approximate analytical methods and computer simulations. 
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APPENDIX A: COLLISIONAL MOMENTS IN THE INELASTIC MAXWELL 

MODEL 

In this Appendix I will derive the collisional moments appearing in the main text. Let 
us consider the general collisional integral 
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A useful identity for an arbitrary function -F(V) is given by 



Irs[F] 



n s Vt d 



w, 



rs 



I 



dvi / rfv 2 / r (V 1 )/ s (V 2 ) J da [F(V'D - F(V0] , 



(A2) 



with 



V'i = Vi - fx sr (l + a rs )(a ■ gi 2 )o\ 



(A3) 
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and gi2 = Vi — V 2 . Now we particularize to F(V) = m r V. Using (A2), one gets 



I rs [m r V] 



n s Vt d 



-m, 



.pi sr (l + a rs ) J dv x J dv 2 f r (V 1 )f s (V 2 ) J der(cr ■ g 12 )er . (A4) 



To perform the angular integration, one needs the result 15 

J da (a- g 12 ) k & = Pk+igi^Sw, 

where 

r(f)r(Mi) 



(A5) 



p k = n d u- l ' r - 



■p / k+d 
" 1 2 



(A6) 



Thus, the integration over a in Eq. (A4) leads to 

I rs [m r V] = -^m r /i sr (l + a rs ) J dv 1 J dv 2 f r (V 1 )f s (V 2 ) (V 1 - V 2 ) 



W r 



p s d 



p sr (l + a rs ) (p s i r - p r j s ) , 



where p r = m r n r and 



j r = J dvm r V/ r (v). 



(A7) 



(A8) 



Next, we particularize to F(V) = m r VV. From the collision rule (A3) it follows that 
yy» _ Vl Vi = p sr (l + a rs ){a • g 12 ) [fi sr (l + a rs )(* • g 12 )<T<5- - (V^ + . (A9) 

To perform the angular integration in (A7), we need now the result 



J da (a- g 12 ) k aa = Y^~ d 9 k i 2 2 (A:gi2gi2 + ff? 2 l) , 



(A10) 



where II is the d x d unit tensor. By using the identity (A2) and Eqs. (A9) and (A10), the 
collisional moment I rs [m r W] can be finally written as 



w r 



/ rs [m r VV] = — -m. 



,p sr (l+a rs ) J dVi J dVa/^VO/.CVa) t^l^lA ( 2gl2gl2 + <&l) 



(Vigi2 + gl2Vi)] 



W r 



p s d 
2 



d + 2 
d 



H sr (l + a rs ) {2p s P r ~ (j r ja + jajr) 

/i sr (l + « rs ) [p s P r + p r P s ~ (jrjs + jsjr) 



+ 



(PrPs + PsPr) ~ jr h 



1 



(AH) 



Here, p r = TrP r /rf = n r T r with 
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J (ivm r VV/ r (v). 



(A12) 



In the simple shear flow state, j r = 0, and the expression (All) for the collisional moment 
I rs [m r W] becomes 



2 w 

J rs [m r VV] = -—-^ul(l + ars f 
d + 2 p s a 



d + 2 \ _ d . 

AW(1 + «„)/ 2 



(A13) 



From Eq. (A13) one can easily get the coefficients A rs and B rs appearing in Eqs. (26)- 
(29). In the case of mechanically equivalent particles (m x = m 2 = m, <j\ = a 2 = cr, an = 
a 22 = a 12 = a), one has Ti = T 2 = T, Pi = P 2 = P, and so the collisional moment 
/ rs [m r VV] = I[mVV] reduces to 



I[mVV] 



1 



-w(l + a) 



[l + a)pl - -(d+ 1 - a)P 



(A14) 



2(d + 2) 

This expression coincides with the one previously derived in the monocomponent case. 19 



APPENDIX B: SOME EXPRESSIONS FOR INELASTIC HARD SPHERES 



The shear flow problem for a binary mixture of IHS was studied in refs. 24 and 25. I give 
here some expressions obtained for this interaction model in the first Sonine approximation 
that have been used along the paper to compare them with the results derived for IMM. 

For IHS the expressions of the coefficients A u , A 2 i, B u , and B\ 2 appearing in Eqs. (25) 
are given by 



y/2 



d(d+2) \a 12 J 



— V 1 #r 1/2 (l + a n ) [a n (d -l) + d+l]v, 



(Bl) 



112 



d(d + 2) 



xiii 12 Q. 



-1/2 ( 9\ + 9 2 



3/2 



(l + ai 2 ) 



29 2 



9 1 + 9, 



+ toi(d- l)(l + ai 2 ) 



v, (B2) 



/^V- 



= 171 , n\ x ± ( " ) H2i0\ l2 {l + an) (2d + 3 - 3a u )u 
d(a + 2) \cr 12 / 



d(d + 2) 



x 2 /i 21 



0i 



1/2 



{[2d - 3/i 2 i(l + a 12 )] (^i + 2 ) + 2(20i + 3# 2 )} i/, (B3) 



d(d + 2) 



X±fJ, 12 



91 



1/2 r 



01 + 02/-, i \ o 

3p 2 i — t (1 + 0:12) - 2 

"2 



(B4) 



Here, Q\ and # 2 are defined by Eq. (21) while v is the effective collision frequency defined 
by Eq. (33). From Eqs. (B1)-(B4) and their counterpart for species 2, the temperature 
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ratio and the elements of the pressure tensor can be easily obtained. In the special case of 
mechanically equivalent particles, Eqs. (B1)-(B4) lead to the following expressions for the 
nonzero elements of the reduced presure tensor P*: 



yy 



p* 

zz 



dd 



d + 1 + (d- l)a 
2d + 3-3a '' 



(B5) 



xy 



A 



d(d + 2)[d + 1 + (d - l)a](l - a) 
2(2d + 3-3a) 2 ' 



(B6) 



(B7) 



Tracer diffusion under shear flow was analyzed in ref. 25 for IHS. The expression of the 
diffusion tensor is 



where the tensor Q, is 



D = P2 (a + ny 1 -Pi, 



(B8) 



(B9) 



Here, A is given by Eq. (56) and the pressure tensor PJj of the excess component is given 
by Eqs. (B5)-(B7). The nonzero elements of P^ can be written in the form given by Eqs. 
(40)-(42) but the functions F, G, and H are 



F = 



d(d + 2) 



3/2 



t*12 



9 3 



1/2 



l + ^I(d-l)(l + 0)(l + a 12 ) 



(BIO) 



1/2 



G 



'^12/^21 



;i + «i 2 ) 



rf(rf + 2) v ^ zl \9(1 + 9)J 
x {2[(d + 2)9 + d + 3] - 3^1(1 + 0)(1 + a 12 )} , 



(Bll) 



H 



3/2 



1/2 



(l + «i 2 ) [3^ 21 (l + ^)(l + « 12 )-2], 



d(d + 2r lz \6(1 + 0)J 
and the reduced shear rate is 

,2 /^2\ 2(d - 1} Pin (l + a 22 )(2d + 3-3a 22 ) 2 (l-a 2 22 ) 



a 12 J d(d+2) 



d+l + (d- l)a 22 



(B12) 



(B13) 
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FIGURES 



FIG. 1. Plot of the reduced elements of the pressure tensor P* as functions of the restitution 
coefficient for a three-dimensional single gas. The solid lines are the results derived here for IMM, 
the dashed lines correspond to the results obtained for IHS from the Sonine approximation, and 
the symbols refer to Monte Carlo simulations for IHS. 24 The dotted line is the result given by Eq. 
(39). 

FIG. 2. Plot of the temperature ratio T1/T2 versus the restitution coefficient in the 
three-dimensional case for n\/n2 = 1, = 1, and two different values of the mass ratio 

m\/m2- m\/rri2 = 2 and m\/mi = 10 . The solid lines are the results derived for IMM, the dashed 
lines correspond to the results obtained for IHS from the Sonine approximation, and the symbols 
refer to Monte Carlo simulations for IHS. 24 

FIG. 3. Plot of the temperature ratio T\jTi as a function of the size ratio o\j(J2 = (mi/ms) 1 ^ 2 
in the two-dimensional case for n\/ri2 = 1 and three different values of the restitution coefficient 
a: a = 0.95, a = 0.9, and a = 0.8. The solid lines are the results derived for IMM, the dashed 
lines correspond to the results obtained for IHS from the Sonine approximation, and the symbols 
refer to Monte Carlo simulations for IHS. 26 

FIG. 4. Plot of the reduced elements of the pressure tensor P* as functions of the restitution 
coefficient in the three-dimensional case for n\/n,2 = 1, 01/02 = E and m\jrwi = 2. The solid lines 
are the results derived for IMM, the dashed lines correspond to the results obtained for IHS from 
the Sonine approximation, and the symbols refer to Monte Carlo simulations for IHS. 24 

FIG. 5. Plot of the reduced elements of the pressure tensor P* as functions of the restitution 
coefficient a in the two-dimensional case for n\/ri2 = 1, and 01/02 = (mi/^) 1 ' 2 = 10. The solid 
lines are the results derived for IMM, the dashed lines correspond to the results obtained for IHS 
from the Sonine approximation, and the symbols refer to Monte Carlo simulations for IHS. 24 

FIG. 6. Dependence of the diagonal and off-diagonal elements of the reduced self-diffusion 
tensor D* on the restitution coefficient a in the three-dimensional case. The solid lines are the 
results derived for IMM while the dashed lines correspond to the results obtained for IHS from the 
Sonine approximation. 

FIG. 7. Dependence of the diagonal and off-diagonal elements of the reduced diffusion tensor 
D* on the restitution coefficient a in the three-dimensional case for o\/o2 = 1 and rxix/m^, = 2. 
The solid lines are the results derived for IMM while the dashed lines correspond to the results 
obtained for IHS from the Sonine approximation. 
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